arXiv:1507.02856vl [astro-ph.GA] 10Jul2015 


Astronomy & Astrophysics manuscript no. multiplydeut © ESO 2015 

July 13, 2015 


Spin-state chemistry of deuterated ammonia 

O. Sipila', J. Harju'’^, P. Caselli*, and S. Schlemmer^ 

* Max-Planck-Institute for Extraterrestrial Physics (MPE), Giessenbachstr. 1, 85748 Garching, Germany 
e-mail: osipila@mpe.mpg.de 

^ Department of Physics, P.O. Box 64, 00014 University of Helsinki, Finland 
^ I. Physikalisches Institut, Universitat zu Kbln, Zulpicher Strasse 77, 50937 Koln, Germany 

Received / Accepted 


ABSTRACT 

Aims. We aim to develop a chemical model that contains a consistent description of spin-state chemistry in reactions involving 
chemical species with multiple deuterons. We apply the model to the specific case of deuterated ammonia, to derive values for the 
various spin-state ratios. 

Methods. We apply symmetry rules in the complete scrambling assumption to calculate branching ratio tables for reactions between 
chemical species that include multiple protons and/or deuterons. New reaction sets for both gas-phase and grain-surface chemistry are 
generated using an automated routine that forms all possible spin-state variants of any given reaction with up to six H/D atoms, using 
the pre-determined branching ratios. Single-point as well as a modified Bonnor-Ebert model are considered to study the density and 
temperature dependence of ammonia and its isotopologs, and the associated spin-state ratios. 

Results. We find that the spin-state ratios of the ammonia isotopologs are, at late times, very different from their statistical values. The 
ratios are rather insensitive to variations in the density, but present strong temperature dependence. We derive high peak values (~ 0.1) 
for the deuterium fraction in ammonia, in agreement with previous (gas-phase) models. The deuterium fractionation is strongest at 
high density, corresponding to a high degree of depletion, and also presents temperature dependence. We find that in the temperature 
range 5 to 20 K, the deuterium fractionation peaks at ~ 15 K while most of the ortho/para (and meta/para for ND 3 ) ratios present a 
minimum at 10 K (ortho/para NH 2 D has instead a maximum at this temperature). 

Conclusions. Owing to the density and temperature dependence found in the abundances and spin-state ratios of ammonia and its 
isotopologs, it is evident that observations of ammonia and its deuterated forms can provide important constraints on the physical 
structure of molecular clouds. 
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1. Introduction 


The chemical evolution of cold molecular cloud cores is char¬ 
acterized by increasing abundances of the nitrogen-containing 
species N 2 H^ and NH 3 , and by the e nhancement of deu terium 
fractionation (e.g., iTine et all IliOOl: iTafalla et all l2002l) . The 
three isotopologs of deuterated ammonia, NH 2 D, NHD 2 , and 
ND 3 , look therefore particularly useful indicators of temporal 
changes in the chemical composition of starless and prestellar 
cores. Besides that the NH 2 D/NH 3 , NHD 2 /NH 2 D, ND 3 /NHD 2 
abundance ratios are like ly to increase with time (iFlower et alJ 
l2006t iRoueff et al.l l2015h . also the ratios of different nuclear 
spin isomers of these molecules can vary. The reactions lead- 
ing to ammonia and its de uterated forms are recapitulated in 
iRodgers & Charni^ (I 2 OOII 1 . Deuteration occurs primarily in 
reactions between NH 3 and deuterated ions. In dense, cold 
clouds this means mainly reactions with the deuterated forms 
of . Through these reactions, the nuclear spin ratios of NH 2 D, 
NHD 2 , and ND 3 become dependent on the nuclear spin ratios 
of H 2 D^, D 2 H^, and DJ. The latter ra tios depend on the H 2 


ortho:para ratio, hereafter 0 /P-H 2 (e.g., IWalmslev etjk 200 ^ 
which also influences the fo rmation rate of ammonia dLe Bourlod 
ll991l:IWirstrom et al.ll 2012 l) . The 0 /P-H 2 r atio is predicted to de- 
crease drastically a s the core evolves (e.g ., IFlower & Wall984 
iPagani et al.l l201ll iBriinken et al.l l2()l4 . While the formation 
rates of different spin isomers of ammonia probably vary be¬ 
cause of competing ion-molecule reactions, grain-surface reac¬ 


tions should produce them in ratios of the nuclear spin statistical 
weights because H and D atoms are added one by one. However, 
the accretion of ammonia molecules synthesized in the gas phase 
leads to gradual changes of the average spin ratios on grains. 


The ground-state rotational transitions of the ortho and para 
species of NH 2 D and NHD 2 , as well as those of the meta and 
para sp ecies of ND 3 are observable with ground-based tele 


scope s dShah & Woottenl200lHFis et al.l2002l ;Jvan_derTaket^ 


Em IRoueff et alJ I2005I: iGerin et all I2006I: Es et al.ri2006l) F 


In the few sources where these have been determined, the 
NH 2 D/NH 3 and NHD 2 /NH 2 D column density ratios have been 
found to lie in the range ~ 0.1 - 0.2, whereas the ND 3 /NHD 2 
is somewhat l ower, ~ 0.02 - 0.1 , implying that ND 3 /NH 3 ~ 
IQ-^ - 10-3 (IRoueff etal. 1 1200.51: iGerin et al.l 12004 ll.is et al.l 
I 2 OO 6 I) . These numbers signify enormous deuterium enhance¬ 
ments in ammonia with respect to the cosmic D/H abundance 
ratio. The observed spin ratios are consistent with the nuclear 
spin statistical weights, taking observational uncertainties into 
account. With the increasing sensitivity of (sub)millimetre tele¬ 
scopes, it will be possible to improve the accuracy of these mea¬ 
surements, and extend the investigation to cover larger and more 
varied samples of dense cores. 


* The fundamental rotation-inversion lines of ortfio-NDs at 615 and 
618 GHz may be observable from the ground in extremely good condi¬ 
tions. 
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A valid description of chemistry in dense cores should be 
able to explain, for example, the observed abundances of the 
different deuterated isotopologs and spin isomers of ammonia, 
which belongs to the most useful probes of dense molecular 
cloud cores. With this goal in mind, we have developed a com¬ 
prehensive gas-grain chemistry model which includes deuterium 
and the spin-state chemistry of molecules containing multiple 
H and/or D nuclei. In an earlier version of this model, the deu¬ 
terium frac tionation and spin selection rules were applied to light 
molecules (ISinila et al.ll2013h . The model has recently been ex¬ 
tended to include heavier deuterated molecules and to take into 
account the spin select ion rules for hydr ogen in the formation 
of ammonia and water (ISinila et al.ll2015l hereafter S15). In the 
present update, we implement the spin selection rules for reac¬ 
tions involving up to six hydrogen and deuterium nuclei. 

As described in Sect. 2 of this paper, the inclusion of the 
spin-state chemistry of multiply deuterated species has required 
a complete revision of the spin separation routine which is now 
based on the symmetries of the nuclear spin functions rather than 
angular momentum algebra used in the previous versions of the 
model. In Sect. 3, we predict the abundances of ammonia and 
its deuterated forms, and the spin ratios of these molecules, as 
functions of time in homogenous models corresponding to dark 
cloud and dense core conditions. In Sect. 4, we discuss the tem¬ 
perature dependence of the deuteration of ammonia, and derive 
radial abundance profiles in two modified Bonnor-Ebert sphere 
models. In this Section, we als o compare our pred ictions with 
the gas phase chemical model of iRoueff et al.l(l2005l) . Finally, in 
Sect. 5 we present our conclusions. An online appendix of this 
paper contains a couple of examples of how to calculate the nu¬ 
clear spin branching ratios using symmetry rules. 


2. Model 


2.1. Theory 

In the present model, the different nuclear spin symmetry states 
of molecules containing two or more hydrogen and/or deuterium 
nuclei are treated as separate species. This implies that for ex¬ 
ample for the doubly deuterated ammonium ion, NH 2 D 2 , we 
keep track of four different forms depending on the two pos¬ 
sible spin symmetries, “ortho” (A) and “para” (B), of the pairs of 
the H and D nuclei. In reactions involving multiple H or D nu¬ 
clei, the branching ratios between different nuclear spin channels 
are assumed to be determined by their pure nuclear spin statisti¬ 
cal weights, i.e., the number of nuclear spin functions belonging 
to each channe l, which are calculated using permutation sym¬ 


metry algebra dOuacklllQTl iPark & LightH2007bl: lHugoll2609t 


iHugo et ^l2009h . The method uses the assumptions of nuclear 
spin conservation and full scrambling of nuclei in reactive colli¬ 
sions. The branching ratios are calculated for reactions involving 

up to six hydrogen or deute rium nuclei. _ 

In our previous papers (ISinila et alJl20I3l : SIS), the nu clear 
spin b ranching ratios were deduced using the method of lOkal 
(I 2 OO 4 I) based on angular momentum algebra. The same method 
w as applied to nucl ear-spin selection rules of nitrogen hydrides 
bv iRist et al.l (l2013l) . For reactions where several deuterium nu¬ 
clei are involved. Oka’s method is, however, impractical, be¬ 
cause species are distinguished by their nuclear spin symmetries, 
but there is no one-to-one correspondence between angular mo¬ 
mentum and sym metry repres entations for systems with three 
or more D nuclei (lHugoll2()(M . A comparison between the to¬ 
tal angular momentum and the total symmetry representations 
of systems containing up to five spin 1/2 and spin 1 particles is 


presented in Table 2.13 of iHugol (l2009h . For multiple H nuclei 
(or any system of spin 1/2 particles), the method using symmetry 
rules gives the same results as Oka’s method, because for these 
systems each value of the spin angular momentum is associated 
with one and only one symmetry representation. Also in the case 
of multiple D nuclei, the total nuclear spin is conserved when us¬ 
ing the symmetry rules. This could be seen from branching ratio 
tables if the symmetry species were subdivided acc ording to the 
total nuclear spin (an example of this is given in iHugol 120091 : 
Table 2.1). 

For the calculation of the statistical branching ratios using 
symmetry rules, one needs to determine the possible symmetry 
species (i.e., representations) of the reactants (F, ( 8 > F^; where ® 
denotes a direct product), the intermediate complex (r„), and the 
products (Fj: ® F/, / -H j = n = k + 1), in appropriate permutation 
groups (Si®S j, Sn, and Sk'^Si). In addition, one needs to derive 
the correlation tables between the group of the complex and its 
subgroups representing the reactants and products. This is done 
applying the standard methods of group theory. The correlation 
tables tell us which nuclear spin symmetry species of the reac¬ 
tion complex are feasible given the symmetries of the reactants, 
and which symmetry s tates of the p roducts are possible when 
the complex decays (cf. lOuacO 19771: Eqs. (12)-(15)). In this de¬ 
scription, the statistical weight of the product F^ ® F; is given 
by /(Fa; ® F;) X dim(rA) x dim(rA), where /(Fa ® F/) is the fre¬ 
quency of Fa ®r; in the total symmetry representation generated 
by the nuclear spin functions of the product molecules in 5 a®> 5/, 
and dim(r) is the dimension of the representation T. Besides be¬ 
ing proportional to the number of nuclear spin functions belong¬ 
ing to each channel, branching ratios calculated in this manner 
are proportional to the number of rovibron ic states occupied by 
these species in the high temperature limit (iPark & Fightj|2007a 
their Appendix A). A couple of examples of the procedure are 
given in the Appendix. The determination of the symm etries of 
the nu clear s pin wave functions is exp lained lucidly in iBunkeil 
dl979h and in iBunker & .lensenl (120061) . In this work, where we 
do not consider the rovibronic states of the molecules, it is suf¬ 
ficient to classify the nuclear spin functions in the complete nu¬ 
clear spin permutation group which does not inclu de the inver¬ 
sion operation, see also ICrabtree & McCalll (l2013l) . Most of the 
statistical branching ratio tables needed in the prese nt study are 
available in Tables III and IV of iHugo et al.l (l2009h . In order to 
cover also reactions with six identical nuclei, like NH 3 -H H 3 and 
ND 3 -H D 3 , and to transfer the spin selection rules to the chem¬ 
ical network, we developed an automated routine for producing 
these tables. The character tables f or permutatio n groups needed 
in this process were obtained using IGA3 (1201 5lPl 


2.2. Application 

We designed a new routine for the spin-state separation of chem¬ 
ical reactions involving multiple protons and/or deuterons, utiliz¬ 
ing the group-theoretical approach described above. As opposed 
to the model described in S15, the new routine produces consis¬ 
tent branching ratios for reactions involving deuterated species, 
instead of copying the spin-state separation rules from the pure- 
hydrogen counterpart reactions. 


^ GAP is a system for computational group theory, 
http://www.gap-system.org 
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Because hydrogen and deuterium are distinguishable, we can 
divide each reaction involving both hydrogen and deuterium into 
two parts. For example, the reaction 


CH 2 D+ H- HDO ^ (COH 3 D+)* 
is divided into 


H 2 DO+ -H CHD 2 


H 2 + H 
D 2 + D 


(H 3 )* 

(D 3 )* 


H 2 + H 

D + D 2 , 


( 1 ) 


( 2 ) 

(3) 


where we have ignored the carbon and oxygen atoms, and the 
charge. The selection rules for hydrogen split reaction (|2l) into 
the following branches: 


oH2-t-H 


pH2-hH 


0 H 2 -t-H 
pH2 + H 
0 H 2 -i-H 
pH 2 + H, 



while the selection rules for deuterium split reaction (O into 

7 

oE )2 + D > D + 0D2 

2 

D -H pD 2 

4 

pE)2 d" D —^ E) -t- oE)2 

5 

D -H pD 2 . 

We then combine the hydrogen and deuterium branches by mul¬ 
tiplying all possible combinations together, and we end up with 
4x4= 16 different branches for reaction O- Note that as a re¬ 
sult, CH 2 D 2 receives a composite spin made up of hydrogen and 
deuterium spin states, for example, C 0 H 2 PD 2 . 

The spin-state separation routine presented here allows us to 
explicitly consider the spin states of all species which contain 
multiple H and/or D atoms. However, we do not track the spin 
states of all species because this would lead to an enormous re¬ 
action network, with a total reaction count far into the hundreds 
of thousands. Instead, we adopt the averaging process described 
in S15, in which we assume a 1:1 spin-state abundance ratio for 
those species whose spin states we do not want to keep track of. 
We have verified that this assumption has negligible influence 
on the results presented below. In this paper, we focus on the 
spin-state chemistry of ammonia and its deuterated forms. We 
also keep the spin states of the species included in the formation 
network of water, so that we can compare our new model to that 
presented in SI5. In reaction ([TJ, the only species for which we 
retain the spin states is H 2 DO^, and the averaging process finally 
results in two branches: 


Fig. 1. Spin-state ratios of NH 3 (upper panels) and H 2 O (lower 
panels) as functions of time at different densities, labeled in the 
upper left panel. The left-hand panels present results calculated 
with the (gas-grain) model of S15, while the right-hand panels 
correspond to the present work. 


states are conserved in the reaction. The reaction -t-H 2 initiat¬ 
ing the am monia synthes i s has b een treated in the same manner 
as done in iDislaire et al.l (1201 2l) and Roueffe^^ (1201 5h . using 
the experimental data from Marquette et akTlOSSh . In addition 
to the updates listed in S15, the new reaction set includes the 
rate coefficients for - 1 - HD and + D 2 from Marquette et 
al., taking into account the slightly different ground-state ener¬ 
gies of 0 D 2 and pD 2 in the reaction enthalpies. For the reac¬ 
tions NH| -H H 2 —> NH | -H H, we adopt the rate coefficients 
from Le Gal et al .) (l2014h . which correspond to those used in 
iRoueff et~ar ( 2005h . 

In what follows, we show results calculated with the gas- 
grain chemical code discussed in S15. In summary, we choose 
the osu_01_2009 reaction set as the basis set which is deuterated 
and spin-state separated according to the approaches discussed 
abov e. The initial grain-sur face chemistry reaction set is adopted 
from ISemenov et an(l 201 (lh . and is similarly deuterated and spin- 
state separated. A complete description of the model including 
adopted model parameters and updates to the (basis) reaction 
sets is presented in S15. 


3. Results 

3.1. Comparison to SI5 


CH 2 D+ -(- HDO 0 H 2 DO+ -H CHD 2 

I 

PH 2 DO+ -H CHD 2 . 

The separation process is fully automated, so any reaction 
involving multiple protons in our base reaction set (see S15) is 
automatically deuterated and spin-state separated according to 
the selection rules discussed above. The only exception to this 
is charge transfer reactions, for which we assume that the spin 


Because of the isomorphism between the symmetry and angu¬ 
lar momentum representations for hydrogenated species, we re¬ 
cover for these the same branch ing ratios in the present model as 
with the method of lOkal (l2004ll employed in S15. However, the 
branching ratios in reactions involving deuterated species are in 
many cases significantly different between the two models. 

In Fig.[T] we compare the o/p ratios of ammonia and water 
predicted by the S15 model and the present model. The tempera¬ 
ture is set to r = 10 K, and the other physical parameters are the 
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same as in S15. The addition of spin states for multiply deuter¬ 
ated species has in general little effect on the spin-state ratios of 
undeuterated species, but some difference is evident, especially 
for water at high density. Inspection of the reaction rates reveals 
that the differences between the models seen in both the ammo¬ 
nia and water o/p ratios are mainly caused by the electron recom¬ 
binations of (o/p) NH 3 D^ and (o/p) H 2 DO^. The S15 model con¬ 
tains some pathways for these reactions that are forbidden by the 
proper selection rules adopted in this work. Nevertheless, the dif¬ 
ference between the models is small, and we have checked that 
this agreement holds for other (undeuterated) species as well. 

The results shown in the following sections were calculated 
with the S15 chemical code, adopting the new reaction sets de¬ 
veloped for this paper. Like in the test above, the temperature is 
in all calculations set to T = 10 K, and the other physical param¬ 
eters are the same as in S15. 

3.2. Deuterium fractionation of ammonia 

By inspecting the production and destruction rates of deuter¬ 
ated forms of ammonia at various stages of the simulation, one 
finds that at late times the deuteration o f ammonia proceeds 
as out lined in the steady-state model of iRodgers & Charnl^ 
(1200ih . The NH 3 abundance has already built up when the gas- 
phase deuteration becomes effective, and deuteron transfer reac¬ 
tions between ammonia and deuterated ions form, e.g., NH 3 D^ 
which can recombine with an electron to give NH 2 D, etc. In 
the present model, however, H 2 D^, D 2 H^, and Dj are the most 
important sources of deuterons at late times. The dominant re¬ 
actions forming NH 2 D, NHD 2 , and ND 3 are shown in Fig. |2l 
This diagram, which is an extension to the schematic NH 3 cy¬ 
cle showmm_Fig^_AJ_^fSf5, should be compared with Fig. 1 
of iRodgers & CharnlevI (1200 ih . Besides the deuteration path¬ 
ways, the diagram shows the constant recycling between NH 3 , 
NHj, and NH^, and similar cyclic processes among their deuter¬ 
ated variants, caused by charge transfer reactions with H^. This 
recycling tends to reduce the deuterium fractionation, because 
NH 2 D^ -h H 2 , for example, leads also to NH^ (we assume a sta¬ 
tistical branching ratios, in this case 4:1 in favour of NH 3 D^). 
The rea ction -H HD —> ND ^ - 1 - H with a very small endother- 
micity dMarauette et al.lflQSSh speeds up the deuteration of am¬ 
monia at early stages, but this pathway has no influence on the 
abundances at the deuterium peak. 

Figure[3] shows the total abundances, i.e., sums over spin 
states, of ammonia and its deuterated forms. At late time at high 
density, ammonia starts to deplete from the gas phase (see also 
S15). For the various forms of deuterated ammonia, an abun¬ 
dance peak is observed which coincides with the initiation of 
deuterium chemistry owing to CO depletion. Like ammonia, its 
isotopologs also deplete from the gas phase, more strongly than 
ammonia itself. This is because of the depletion of deuterated 
molecules (including the ammonia isotopologs) onto grain sur¬ 
faces which traps deuterium onto the surfaces; this results ulti¬ 
mately in HD depletion and con sequently the over all gas-phase 
deuteration efficiency decreases (ISinila et al.ll2013l : S15). 

The abundance ratios of ammonia and its isotopologs are 
presented in Fig.H] The deuterium fractionation is highest at 
high density where depletion is strongest. The relative depletion 
of the ammonia isotolopogs is rather uniform, with all abun¬ 
dance ratios ultimately dropping by a factor of 3-4 from their 
peak values. A high degree of fractionation is observed, and the 
peak values for all of the ratios are > 0.1. This issue is further 
discussed in Sect. 14.31 where we com pare our gas-grain m odel- 
ing results to the gas-phase models of iRoueff et alJ([^05h . 


Table 1. Spin symmetries of molecules with 2-4 H or D nu¬ 
clei. The numbers in front of the symmetry labels indicate their 
frequencies in the total nuclear spin representation of the com¬ 
plex. The values of the possible nuclear spin angular momenta, /, 
and the nuclear spin statistical weights, gj, are given below each 
symmetry species, as well as their ’’para”, ’’ortho”, and ”meta” 
appellations, when applicable. 


Hz 

3A B D, 

6A 3B 

H 3 

/ = 1 / = 0 
ft = 3 ft = 1 
ortho para 

4Ai 2E Di 

7 = 0,2 7=1 

g; = 6 g, = 3 
ortho para 

lOAi Az 8E 

H 4 

7 = 3/2 7=1/2 

ft = 4 g; = 4 

ortho para 

5Ai E 3fi 

7=1,3 7 = 0 7=1,2 

g/ = 10 g, = 1 g; = 16 

meta para ortho 

D 4 

1=2 7=0 7=1 
g/ = 5 g/ = 2 g/ = 9 
meta para ortho 

15Ai 6E 15Fi 3 E 2 


7 = 0,2,4 7 = 0,2 7 = 1,2,3 7 =1 
ft = 15 g/ = 12 g; = 45 g; = 9 




t w 


Fig. 3. Abundances of ammonia and its deuterated forms as func¬ 
tions of time at different densities, labeled in the upper right 
panel. 


3.3. The spin-state ratios of deuterated ammonia 

In Fig.|5] we show the spin-state ratios of the deuterated forms 
of ammonia. All of the ratios show dependence on both time 
and density. For all of the ammonia isotopologs, the most im¬ 
portant formation mechanism is the electron recombination of 
the (singly or multiply) deuterated ammonium ion. To comple¬ 
ment Fig.|5j we show in Fig.the spin-state ratios of selected 
ions. In the following, we quote the spin symmetries of the vari¬ 
ous species; these are tabulated in Table[T] 

At late times, the ortho:para ratio of NH 2 D is largely deter¬ 
mined by the cycle NH 2 D ^NH 2 D+^NH 3 D+^NH 2 D (see 
Fig. 12). The spin ratio of NH 3 was discussed in Appendix A 
of S15. The electron recombination of NH 3 D^ yields ortho and 
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Fig. 2. Principal reactions determining the abundances and spin ratios of deuterated ammonia isotopologs at late stages of chemical 
evolution in cold, dense cores. The reactions NH 3 _xDx + H 3 _yDy are indicated with blue arrows. 



Fig. 4. Abundance ratios of ammonia and its deuterated forms as functions of time at different densities, labeled in the left-hand 
panel. 


para NH 2 D with the following branching ratios: 

oNH 3 D+-i-e~ 0 NH 2 D-HH 

, '-h 

PNH 3 D+ -h e ^ PNH 2 D -h H 

'-h 

^ 0 NH 2 D-HH. 

Since the destruction rates of 0 NH 2 D and PNH 2 D (mainly 
through reactions with and Hp are similar, we get o/p- 
NH 2 D ~ 1 -H 2 o/p-NH 3 D^. Moreover, because the o/p ratio of 
NH 3 D^ is ~ 0.4 - 0.5 at late times depending on the density 
(Fig.| 6 ]l, 0 /PNH 2 D should obtain the values 2 (low density) or 
1.8 (high density), in agreement with Fig.|5] 

The o/p ratios of NH 2 D and NH 3 D^ are, however, closely 
coupled because many of the most important destruction path¬ 
ways of NH 2 D return NH 3 D^. This can be seen by inspec¬ 
tion of Figs.|5] and which show the results of the full simu¬ 
lation, and from the schematic presentation of the situation at 
late times in Fig. |2] The spin state is conserved in the charge 

H+ 

exhange reaction NH 2 D —»NH 2 D^. Assuming total scrambling 
of nuclei in the reaction NH 2 D^ -H H 2 , oNH 2 D^ -H pH 2 gives 
rise to the intermediate complex NH 4 (Fi)D^ which dissociates 
to oNH 3 D^ and pNH 3 D^ in ratio 1:2 (the symmetry of H 4 is 


indicated in the brackets, see Table III of iHugo et al.l 1200^ . 
The reaction pNH 2 D^ -H pH 2 produces only pNH 3 D^ through 
the intermediate complex NH 4 (£’)D^ (only reaction with pH 2 
are considered because the 0 H 2 abundance is very low at this 
stage). Using these branching ratios one obtains that the closed 

cycle NH 2 D ^ NH 2 D+ ^ NH 3 D+ NH 2 D should result in o/p- 
NH 2 D ~ 2.3. As discussed above, in the full reaction network 
this ratio is, however, about 20 % lower. 

The doubly and triply deuterated froms of ammonia are pri¬ 
marily formed from the electron recombination of NH 2 D 2 and 
NHD 3 , respectively. (The recombination of the fully deuter¬ 
ated ammonium ion ND^ contributes about 10% to ND 3 .) The 

cycle NHD 2 ^ NH 2 D^ S NH 2 D 2 ^ NHD 2 preserves the spin 
state of D 2 . The same is true for D 3 in the corresponding cycle of 
ND 3 . NH 2 D 2 forms mainly from the reaction D 2 H^ + NH 3 —> 
NH 2 D 2 + H 2 , where the spin state of D 2 H^ is conserved. NHDj 
receives contribution from both D 2 H^ and through the reac¬ 
tions D+ H-NH 3 ^ NHD+ -t-H 2 and D 2 H+ -hNH 2 D ^ NHD+ H- 
H 2 . In the first of these reactions, the spin ratios of are just 
copied to ND 3 . For the s econd reaction, the branching ratios can 
be read from Table IV of iHugo et al.l(l2009h . For example, taking 
the most common species, the reaction oD 2 H^ + 0 NHD 2 forms 
an intermediate reaction complex (NH 2 DP* with Ai, E, and Fi 
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Fig. 6 . Spin-state ratios of the deuterated forms of the ammonium ion, as well as D 2 H^ and D 3 , as functions of time at different 
densities, labeled in the lower left panel. 




Fig. 5. Spin-state ratios of the deuterated forms of ammonia as 
functions of time at different densities, labeled in the lower left 
panel. 


symmetries of D 4 in statistical ratios 15:6:15, and, taking these 
ratios into account, the complex decay into NHDJ + H produces 
A I (“meta”) and E (“ortho”) symmetries of D 3 with probabilities 


5/9 and 4/9, respectively (the A 2 “para” symmetry is not formed 
in this reaction because in D 4 ^ D 3 -f D, A| Ai, E —> E, and 
El ^Ai+E, Table IV of iHugo et alJl200l . 

What is said above implies that the spin ratios of NHD 2 and 
ND 3 are tied to those of D 2 H^ and , and this is also evident 
from Figs. |5] and | 6 ] The o/p ratio of NHD 2 increases constantly 
until very late times, levelling off only after 10 ^ years in the mod¬ 
els shown in Fig.|5](i.e., outside the range shown here). This de¬ 
pends on the D 2 H^ o/p ratio which increases strongly as a func¬ 
tion of time (Fig.| 6 ]l. From Fig.|5]it is observed that para-ND 3 is 
clearly the least abundant of the three spin states of ND 3 , while 
the ortho and meta states are equally abundant at late times re¬ 
gardless of the density. This behaviour reflects the spin ratios of 
Dj shown in Fig.| 6 l 

The statistical values for the spin-state ratios of ortho/para 
NFI 2 D, ortho/para NFID 2 and ortho/meta/para ND 3 are 3:1, 2:1, 
and 16:10:1, respectively. From Fig.|5]it is evident that the var¬ 
ious ratios deviate strongly from the statistical values at late 
times, as was already found for ammonia in S15 (see also Fig.ITJ. 
The various spin-state ratios are generally rather insensitive to 
variations in the density. Some density dependence is usually 
present between 10 ^ and 10 ® years of chemical evolution, but 
even in this time interval the difference between the low and high 
density cases is less than a factor of 2 . 

Assuming nuclear spin conservation and that the same selec¬ 
tion rules as used in the gas phase can be applied to the H and 
D addition reactions on grain surfaces, the different spin isomers 
of ammonia should form on grains according to their nuclear 
spin statistical weights indicated above. This is what is seen in 
the beginning of the simulation. Later on, when the gas-phase 
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ammonia starts to accrete, the spin ratios change also on grains, 
finally mimicking the ratios in the gas. 

Finally, we note that for the lowest density considered (tin - 
10 ^ cm“^), many of the plotted spin ratios are increasing even at 
very long timescales, for example the m/p ratio of ND 3 which 
increases by a factor of ~ 2 from 10 ® to 10 ^ yr at this density. 
We note that this effect is due to the very low density and the 
spin ratios are expected to stabilize at still longer timescales. 

3.4. Temperature dependence of ammonia D fractionation 
and spin-state ratios 

It is well established that deuterium fractionation depends not 
only on the density, but also on the temperature. In particular, 
deuteration is expected (in starless cores) to be at its strongest in 
the temperature range ~ 10 - 20 K, because on the one hand this 
range lies below the sublimation temperature of CO, allowing 
the main driver ion of deuterium chemistry, H 2 D^, to transfer its 
deuteron to other species, and on the other hand deuteration effi¬ 
ciency drops toward low temperatures because of the decreasing 
accretion rates (see below). 

We plot in Fig.Qthe abundances of ammonia and its deuter¬ 
ated forms, and the associated abundance ratios, at different 
temperatures. Evidently, the deuterium fractionation is strongest 
around T = 15 K; both the abundances and the deuterium frac¬ 
tionation peak at this temperature. An analysis of the reaction 
rates shows that if the temperature is reduced from 15K to 
5 K, the decreased accretion rates lead to relatively high abun¬ 
dances of atomic H and D in the gas phase, which enable effi¬ 
cient conversion of D 2 H^ to H 2 D^ through the exchange reac¬ 
tion D 2 H^ -h H —» H 2 D^ -t- D. This decreases the ammonia deu¬ 
terium fractionation because D 2 H^ is an important part of the 
production paths of the deuterated forms of ammonia (Fig.|2]i. 
On the other hand, as the temperature is increased from 15 K 
to 20 K, CO starts to desorb efficiently, hindering deuteration in 
general. 

The spin-state ratios of the ammonia isotopologs at different 
temperatures are shown in Fig.| 8 ] It is observed that the ratios 
change significantly as the temperature is varied, and that the 
variation is much larger than that induced by changes in the den¬ 
sity. A minimum in the ortho/para (and meta/para for ND 3 ) ratios 
is observed at T = 10 K, except for the ortho/para ratio of NH 2 D 
which has a maximum at this temperature. These are in contrast 
to the overall deuteration, which peaks at T = 15 K. 

As explained in Sect. l3.3l the spin-state ratios of ND 3 are de¬ 
termined by those of NHD^, which are in turn determined by 
D 2 H^ and D 3 . Moving from 10 K to 20 K, the ortho/para ratio of 
H 2 increases, leading to higher ortho/para D 2 H^ and meta/ortho 
(and meta/para) DJ ratios, which ultimately increase the or¬ 
tho/para and meta/para ratios of ND 3 with respect to the values at 
10 K. On the other hand at 5 K, the H 2 ortho/para ratio is low, but 
the abundances of atomic H and D increase owing to decreased 
accretion rates, and the proton/deuteron exchange reactions with 
the H 3 isotopologs (for example, D 2 H^ -t- D —> DJ - 1 - H) in¬ 
crease in importance in determining the various spin-state ratios. 
The spin-state ratios of NH 2 D and NHD 2 are affected by similar 
effects with varying temperature. 

4. Discussion 

4.1. Radial profiles for deuterated ammonia 

From the above it is clear that both the deuterium fractiona¬ 
tion of ammonia and the spin-state ratios depend on time and 



Fig. 8. Spin-state ratios of the deuterated forms of ammonia as 
functions of time at different temperatures, labeled in the lower 
left panel. A constant density of uh = 10® cm ® is adopted. 


on the physical parameters (density, temperature). To investigate 
the implication to observations of deuterated ammonia, we con¬ 
structed radial abundance profiles for ammonia a nd its deuter¬ 
ate d forms by followi ng the method discussed in ISinilal (12012h 
an d lSinila etal.l(l2()H . in which chemical calculations are com¬ 
bined with a radiative transfer model to produce self-consistent 
chemical abundance profiles as functions of density, temperature 
and time. For these calculations , we adopt as the core model a 
modi f ied Bonnor-Ebert sphere (iGalli et al.l l2002t iKeto & Eieldl 
I 2 OO 5 I: [^pila et al.ll201 Ih of mass M — 1.0 Mq. We assumed a 
visual extinction Ay = 10 mag at the edge of the core. 

The density and temperature profiles of the model core are 
shown in the left-hand panel in Eig.|9] The density and dust tem¬ 
perature profiles hardly change as functions of time. The density 
is not high enough even in the core center to equalize the dust 
and gas temperatures. On the other hand, the gas temperature de¬ 
creases toward the outer areas because of efficient line cooling 
by (mainly) CO. As the chemical evolution progresses, CO de¬ 
pletes on grain surfaces, reducing the cooling rates, and the gas 
temperature rises. 

The calculated radial abundance profiles are shown in the 
middle and right-hand panels in Eig.|9] The time steps shown 
in the figure refer to years of chemical evolution with respect 
to an initially atomic initial state (except for H and D which 
are initially in H 2 and HD; see S15). The abundances of ammo¬ 
nia and its deuterated isotopologs generally increase as a func¬ 
tion of time. Notably, the abundance of ND 3 increases by about 
three orders of magnitude across the core from f = 5 x 10"^ yr to 
t - 1 X 10 ® yr, and the increase in time is stronger for species with 
more D atoms. The difference in abundance between the core 
center and edge can be close to two orders of magnitude depend¬ 
ing on the species considered and on the time. Depletion slightly 
decreases the abundances at the core center at the last time step, 
while the abundances keep increasing at the edge where the den¬ 
sity is lower and little depletion is present. An exception to the 
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Fig. 7. Upper panels: abundances of ammonia and its deuterated forms as functions of time at different temperatures, labeled in the 
left-hand panel. Lower panels: abundance ratios of ammonia and its deuterated forms as functions of time at different temperatures. 
In all panels, a constant density of «h = 10^ cm“^ is adopted. 



Fig. 9. Left-hand panel: Gas density (black), gas temperature (red) and dust temperature (blue) as functions of distance from the 
core center at different time steps (indicated in the middle panel). Middle panel: Abundances of NH 3 and NH 2 D as functions of 
distance from the core center at different time steps. Right-hand panel: Abundances of NHD 2 and ND 3 as functions of distance from 
the core center at different time steps (same legend as in the other panels). The panels correspond to an MBES of mass 1.0 Mq. 


latter behavior is ammonia, whose abundance has a minimum in 
the outer core around f = 5 x 10^ yr; this kind of behavior is also 
displayed for example in Fig.|2]and is discussed in Sect. l4.2l 

The results shown here demonstrate on the one hand the im¬ 
portance of accurate knowledge of the physical structure, and on 
the other hand the necessity of detailed chemical modeling in as¬ 


certaining the chemical age of the object in question. However, 
quantifying the effect of the varying abundance profiles on ob¬ 
servable emission requires calculation of the line profiles with a 
radiative transfer model; this will the subject of a future paper. 
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Fig. 10. Fractional abundances (with respect to hh) of selected 
species (labeled in the plot) in a homogeneous model with iiu - 
lO^cm ^ r = lOK. 


4.2. Evolution of the ammonia abundance 

Characteristic for the time dependence of the NH 3 abundance at 
high densities, and to a lesser extent for that of NH 2 D, is a mini¬ 
mum occurring somewhere between 10 ^ and 10 ® yr, and a rapid 
recovery after this minimum. This tendency is most clearly vis¬ 
ible in Fig.|7] The behavior is very different from that observed 
in the case of CO or CS, for example, but resembles the varia- 
tion of the abundan ces of HCN and HNC seen in the models of 
iLoison et al.l (l2014ll . 

The ammonia production is assumed to start slowly, the 
widely recognized bottle-neck being the slightly endothermic re¬ 
action + H 2 —> NH^ -f H which in practice works only with 
0 H 2 in cold clouds (e.g., Dislaire et al. 2012; Roueff et al. 2015). 
Once NH^ is formed it is quickly converted to ammonia through 
NH^ + H 2 —» NHj -h H etc. In fact, the flow through the bottle¬ 
neck is so slow that most of the time the chemical network may 
find other means to produce NH^. In the present model, NH^ is 
copiously formed at early times (~ 10 "^ yr) from the dissocia¬ 
tive ionization of HNC by He^, the bottle-neck reaction coming 
in second. The efficiency of this channel decreases along with 
the depletion of HNC, and because NH 3 accretes onto grains as 
well, its gas-phase abundance decreases some time before 10 ® 
yr, or slightly after that depending on the density. 

A strong increase in the ammonia abundance is seen, how¬ 
ever, after the minimum, as an indirect consequence of the de¬ 
pletion of CO and neutrals like HCN and OH which are ion¬ 
ized in charge transfer reactions with H^. The disappearance 
of these molecules enhances greatly the abundances of H^ and 
H^. Firstly, Hj converts part of N 2 to NH through the reaction 
H;-hN 2 ^ N 2 H+-(-H 2 , followed by N 2 H+ -He- ^ N 2 -hH/NH-)- 
N. Although less than 10% of the recombinations lead to NH, 
this channel is important for ammonia, because the increased 
abundance of H^ directs NH to the ammonia track through the 
charge transfer reaction H^ -H NH —> NH^ + H. Likewise, all 
the recombination products of NH 3 + and NH^ + (NH, 
NH 2 , NH 3 ) are efficiently returned to the ammonia cycle by H^. 
In the present model, ammonia passes several times the cycle 



Fig. 11. Abundance ratios of ammonia and its deuterated forms 
(labeled in the plot) as functions of time in a gas-phase chemical 
model, for different initial values of the H 2 ortho/para ratio: 3 
(black lines), IQ-® (blue), or IQ-® (green). 

NH 3 ^ NHj —> NH 4 ^ NH 3 described in Sect. l3.2l before col¬ 
liding with a dust grain. In Fig.[T0]we show the fractional abun¬ 
dances of NH 3 , NH, H^, Hj, and HNC as functions of time for 
a homogeneous model with tin = 10® cm ®, T = lOK. In this 
model, NH 3 and NH are closely correlated after 10® yr and fol¬ 
low the gradual decrease of the N 2 abundance (not shown) at 
very late times. 

4.3. Comparison to lRoueff et al\ f2005l) 

iRoueff et al.l (l2005h presented observations of NH 2 D, NHD 2 , 
and ND 3 toward prestellar cores, and discussed the interpretation 
of the observations in the context of a gas-phase chemical model. 
Their model is to our knowledge the most directly comparable 
to our own, in terms of the chemistry of deuterated ammonia. 

Our model includes time-dependent depletioifl and a chem¬ 
ical steady-state is not reached in the time intervals consid- 
ered here. However , comparing the values given in Table 9 in 
iRoueff et al.l (l2005h with our results near the deuterium peak 
for each density (Fig.O, we find that the agreement between 
the models is very good, despite the fact that we are using a 
gas-grain model with different initial abundances and a different 
value for the cosmic ray ionization rate (here 1.3 x 10 '® s ®; 
2 X 10 '® s ' in Roueff et al.). It should be noted, however, that 
in our model, depletion occurs also at a density of 10 ^ cm ®, 
whereas in the Roueff et al. model no depletion is assumed at 
this density. 

As an additional check, we ran another calculation using a 
gas-phase model (i.e., no depletion) at «(H 2 ) = 10 “® cm-®, for 
three different values of the initial H 2 ortho/para ratio. The re¬ 
sulting abundance ratios for ammonia and its deuterated forms 
are shown in Fig.[TT] The gas-phase model predicts a lower deu¬ 
terium fractionation than the gas-grain model, even for very low 
initial H 2 ortho/para ratios. The main reason for this is that the 

^ In o ur model, nitrogen a nd deuterium deplete as well, unlike in the 
model of lRoueff et all l l2005ll . 
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high abundance of CO hinders the formation of H 2 D^, D 2 H^ and 
D 3 , which are in our model the principal agents in the formation 
of deuterated ammonia. 

iRoueff et alJ (l2005h assumed that the release of hydrogen is 
favored (in a 2:1 ratio) over the release of deuterium in the elec¬ 
tron recombinations of the various deuterated ammonium ions, 
whereas in this work the H/D branching ratios for the recombi¬ 
nation reactions are statistical. We tested the effect of assuming 
a 2:1 ratio for the recombination reactions and found that this 
has only a marginal effect on the abundances of ammonia and its 


isotopologs. We n ote that we adopt the NH^ + H 2 


NH+ -i-H 


reaction rates f romiLe Gal et al. ( 20141) whi ch correspond to th e 
values used by iRoueff et~ar ( 20051: see also Roueffet^ 2015h . 
We tested the new coefficients given by IRoueff et al. ( 2015h for 
these reactions and found that they have only a minor effect on 
our results. 

The differences between our gas-phase model and that of 
IRoueff et al.l (l2005h can perhaps be explained by the adoption of 
different rate coefficients for the - 1 - H 2 reacting system, since 
the contribution of the isotopologs to the chemistry of deuter¬ 
ated ammonia is significant as discussed earlier in this paper. 
Finally we reiterate that the depletion present in our gas-grain 
model even at lower densities enhances d euteration and bring s 
the results close to the gas-phase model of IRoueff et al.l (l2005lf 


5. Conclusions 

We calculated branching ratios for reactions involving multiply 
deuterated species using group theory. This method produces 
identical branching ratios for reactions without d euterium, w hich 
can be derived with angular momentum algebra (I0kal2004l) . Our 
new model allows us to track the spin states of any species that 
includes multiple protons and/or deuterons. However, in this pa¬ 
per, we concentrated on the spin-state chemistry of deuterated 
ammonia. 

We found that the deuterium fractionation of ammonia is 
high especially at high density: the NH 2 D/NH 3 , NHD 2 /NH 2 D, 
and ND 3 /NHD 2 all peak at ~ 10% for «h = lO^cm^^. The 
hi gh abundance ratios are in line with earlier gas-phase models 
of IRoueff et al.l (l2005l) . The various spin-state ratios were found 
to deviate strongly from the statistical values especially at late 
times, but on the other hand the ratios appear rather insensitive 
to variations in the density. We also found that the abundances of 
ammonia and its isotopologs and the associated spin-state ratios 
are highly dependent on the temperature; the meta/para ratio of 
ND 3 for example can vary by a factor of ~ 5 in the temperature 
range 5 - 20 K. 

The abundances of NH 2 D, NHD 2 , and ND 3 , and their spin- 
state ratios in the gas phase, are determined by competing ion- 
molecular reactions, the rates of which depend strongly on the 
kinetic temperature. The deuterium fractionation is most effi¬ 
cient at 15 K. The temperature dependences of the spin-state ra¬ 
tios are less straightforward: the para species of NHD 2 and ND 3 
reach their abundance maxima at 10 K, wheras para-NH 2 D has 
a minimum at this temperature. The reactions leading to deuter¬ 
ated ammonia involve H 2 D^, D 2 H^, and , and benefit from 
the disappearance of CO which occurs faster at high densities. 
Therefore, for objects with density/temperature gradients, such 
as starless and prestellar cores, significant variations as func¬ 
tions of radius and time can be expected in the abundances and 
the spin-state ratios of NH 2 D, NHD 2 , and ND 3 . Sensitive obser¬ 
vations of these molecules, combined with chemical and radia¬ 
tive transfer modelling can provide stringent constraints on the 


physical and chemical structures, and the evolutionary stages of 
molecular clouds. 

The present model does not include corrections to rate co¬ 
efficients resulting from differences in the zero-point energies 
of deuterated and undeute rated specie s , exce pt for the Hj -t- H 2 
reacting system given by Hugo et ^ (l2009l) and the + H 2 
reaction (iMarauette et al.lll988l) . MOTeover. the use of pure nu¬ 
clear spin statistical weights is based on the implicit assumption 
that the exothermicity of the reaction is sufficiently large so that 
Pauli allowed rotational levels of each product species can be 
populated. For the principal reactions contributing to the inter- 
stellar ammonia , this is likely to be true (see also discussion in 
iRist et alll2013h . However, it is clear that the present model is 
still incomplete, and an in-depth study of energetics is warranted. 
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Table A.l. Statistical nuclear spin branching ratios for reactions 
H 2 -f H ^ H 3 and D 2 -H D ^ D 3 . 


H 3 

Hz 

A (8 A 

-rH 

B^A 

2 

Da 

Dz 
A 8A 
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4Ai 

4 

0 

4 

lOAi 

10 

0 

10 
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0 

0 

0 

IA 2 

0 

1 

1 

2E 

2 

2 

4 

8B 

8 

8 

16 

z 

6 

2 

8 

2 

18 

9 

27 


Appendix A: Calculation of nuclear spin branching 
ratios using symmetry rules 

The method used in the derivation of the nuclear spin branch¬ 
ing ratios is demonstrated by three examples, all dealing with 
varieties of the reaction H++ NH 3 ^ (N H+)* NH| -H H 2 . 
The me thod is described in lOuackf (Il97^ and in iPark & Lighll 
(l2007bl) . The collision partners are classified, or labeled, accord¬ 
ing to the irreducible representations of an appropriate permuta¬ 
tion group. These labels tell how the symmetrized nuclear spin 
functions belonging to this species transform under operations of 
the group, like the transposition of tw o identical nuclei and th e 
permutation of n identical nuclei (e.g.- lBunker & .Tensenll200^ . 
The branching ratios are obtained from correlations between 
the group of the intermediate reaction complex and th e direct 
produ ct groups representing the reactants and products (lOuackI 
Il977h . The statistical weight of each channel corresponds to the 
number of nuclear spin functions belonging to that symmetry 
species. 

A. 1. Three identical nuclei 

As the first example we consider the reaction D 2 H^ - 1 - NH 2 D —» 
(NH 3 DP* —» NHDj + H 2 . The H and D symmetries are treated 
separately. The same result would be obtained by considering 
the mixed system in appropriate direct product groups. From the 
hydrogen point of view, H 2 and H form H 3 which dissociates 
back to H 2 -t- H. On the deuterium side, the reaction forms D 3 
from D 2 + D. In the present model we assume that the branching 
ratios are equal to the pure nuclear spin symmetry induction and 
subduction statistical weights, which for these r eactions can be 
read from Tables III and IV of iHugo et al.l (l2009h . We reproduce 
the appropriate tables in Table lA.ll 

The numbers in the leftmost columns of Table lA.ll are the 
frequencies of the irreducible representations Ai, A 2 , and E, in 
the representations of the permutation group S 3 generated by 
the nuclear spin functions of and . This means that the 
symmetry representations of H 3 and D 3 in 5 3 can be reduced 
into the following direct sums (denoted by ©): 

r"-*(H3) = 4Ai©2£ , r‘’'^(D3 ) = 10Ai© 1A2®8£ (A.l) 

(e.g., iBunker & JensenI |2006t iHugo et ^l2009l) . According to 
Mulliken’s notation used here, one-dimensional representations 
are labelled either A or B, two dimensional are labelled E, three- 
dimensional F, four-dimensional G, etc. The total nuclear spin 
statistical weight of each symmetry species can be read from 
these decompositions by multiplying the frequency by the di¬ 
mension of the representation, /(T,) x dim(r,). The Greek ap¬ 
pellations, para, meta, and ortho, are used to indicate these sta¬ 
tistical weights for the simplest molecules, so that “ortho” has 
the highest statistical weight. The “para” (E, I - I /2) and “or¬ 
tho” (Ai, / = 3/2) species of have, however, equal statistical 
weights. has no nuclear spin functions of species A 2 which 


is antisymmetric with respect to interchange of two nuclefl. H 2 
and D 2 can have the symmetries A or B, and the single nuclei H 
and D have the A symmetry. The number in the bottom right cor¬ 
ner of each table gives the total number of linearly independent 
nuclear spin functions of the molecule, (2/ + 1)^, where I -Xjl 
for H 3 and 7=1 for D 3 . The branching ratios in Table lA.fi are 
obtained by multiplying the correlation table between S 3 and the 
product group 52 ® 5 1 by appropriate frequencies. This correla¬ 
tion table is given by 


S 3 

52 ® 5i 
A(g)A B®A 

Ai 

1 

0 

A 2 

0 

1 

E 

1 

1 


The correlation table 5 3 <-> 5 2 ® 5 1 presented above is de¬ 
rived using the standard methods of group theory, i.e., by form¬ 
ing the character table of direct product group ^2 ® 5 1 , and the 
character table for the subduced representations of 5 3 on 5 2 ® 51 , 
S 3 i 5 2 ® 5 1 . The latter is obtained from that of 5 3 by strip¬ 
ping it from operations which are unfeasible in the direct prod¬ 
uct group 52 ® 5 1 , i.e., the class of cyclic permutations of three 
nuclei, {(123), (132)}. Finally, the orthogonality relation between 
the characters of irreducible representations is used. 

The Tables IATTI can be read either from left to right (’’sub¬ 
duction”) or top down (’’induction”). For example, one can see 
that the Ai species of H 3 and D 3 dissociate exclusively to the A 
symmetry of H 2 or D 2 , whereas the E species dissociate to both 
A and B symmetries (“ortho” and“para”, respectively) of H 2 or 
D 2 with equal probability. In the reaction 0 D 2 H- D (A (g) A), the 
Ai (“meta”) and E (“ortho”) species of D 3 are formed with the 
probabilities 10/18 and 8/18, respectively, whereas in the reac¬ 
tion pD 2 + D (B ® A), the A 2 (“para”) and E (“ortho”) species of 
D 3 are formed with the probabilities 1/9 and 8/9. These branch¬ 
ing ratios determine relative abundances of NH 0 D 3 , NHmD^, 
and NHpDj resulting from the reaction D 2 H^ -H NH 2 D. 

A.2. Six identical nuclei 

In order to determine the nuclear spin branching ratios of the re¬ 
actions H+ -t- NH 3 ^ (NH+)* ^ NH| -H H 2 and D+ H- ND 3 ^ 
(NDg)* —» ND^ -H D 2 one needs to determine the nuclear spin 
symmetry species of He and De in the permutation group 5 e. We 
calculate first the characters, ;if(7?), of the representations gener¬ 
ated by the nucle ar spin functions of He an d De under the op¬ 
erations R of 5e. ICrabtree & McCalil (1201 3l) give a convenient 
formula for this: 

;^(7;) = (2/+l)"-', (A.2) 

where I is the nuclear spin (1/2 for H and 1 for D), 11 is the 
number of identical nuclei (here 6 ), and x is the number of trans¬ 
positions or pairwise permutations which correspond to the op¬ 
eration R. The number of interchanges corresponding to a cyclic 
permutation of k elements is k -1, and so x - k-l. For example, 
X - 2 for the ternary permut ation (123). T he character table of 
56 is easily obtainable using Iga 3 (l2015l) . By applying the or¬ 
thogonality relation of the characters of irreducible representa¬ 
tions, we derive the following decompositions for the symmetry 
representations of He and De: 

F"'" (H 6 ) =7Ai®5 77i®17/4®3Li , (A.3) 

The fulfilment of the symmetrization postulate for each nuclear spin 
species is taken care of by the complete internal wavefunction consist¬ 
ing of the nuclear spin, ro-vibrational and electronic parts. 
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Table A.2. Statistical nuclear spin branching ratios of the reac¬ 
tion H6 ^ H3 -I- H3. 


Hfi 

Ai, Ai 

H3 -hHs 
A,,F F,A, 

F,F 

2 

7Ai 

7 

0 

0 

0 

7 

5F, 

5 

10 

10 

0 

25 

IH 4 

1 

0 

0 

4 

5 

3Li 

3 

6 

6 

12 

27 

2 

16 

16 

16 

16 

64 


Table A.3. Statistical nuclear spin branching ratios of the reac¬ 
tion He ^ H 4 -H H 2 . 


Hfi 

Ai,A 

Ai,B 

H4 

F,A 

+ H2 

E,B 

Fi,A 

Fi,B 

2 

7Ai 

7 

0 

0 

0 

0 

0 

7 

5 Hi 

5 

5 

0 

0 

15 

0 

25 

IH 4 

0 

0 

0 

2 

3 

0 

5 

3Li 

3 

0 

6 

0 

9 

9 

27 

2 

15 

5 

6 

2 

27 

9 

64 


-r NH3 (NH+)* NH; -r H2 



^"*•(06) = 28Ai®35//i®1//3®10//4®27Li®10Mi®85 .(A.4) 

The correlation tables ^ Sj,®Si and S(, ^ 84'^ S2 needed 
for the induction and subduction statistics are derived as outlined 
in the previous example. Finally, the branching ratio tables listed 
in Tables lAT2llA.31lA~4l and lA.SI are obtained by multiplying the 
correlation tables by the frequencies of irreducible representa¬ 
tions appearing in the decompositions given above. In addition, 
the columns of the correlation table are multiplied by the di¬ 
mensions of the subduced representations in order to make the 
weights correspond to the total number of nuclear spin functions 
belonging to each representation. In the previous example, all 
products had one-dimensional representations. The correct nor¬ 
malisation of the weights can be checked by inspecting the last 
rows and columns of the tables. In these tables, species with zero 
frequencies are not listed. 

The use of the branching ratio tables is illustrated in 
Figs. IA.II and IA.21 Fig. lA. II gives the outcome of the reaction 
Hj + NH 3 assuming that one of the molecules is of the Ai (“or¬ 
tho”) symmetry and the other one is E (“para”). The branch¬ 
ing probabilities indicated next to the arrows are obtained from 
Tables IAT 2 I (reading it top down) and from Table |A3] (reading it 
from left to right). After combining the probabilities of pathways 
leading to different symmetry species of NH^ one finds that the 
F\,A\, and E species are formed in ratios 15:7 : 2 in this reac¬ 
tion. Fig. lA.2l describes the branching of the reaction Dj - 1 - ND 3 
assuming that both species are Ai (“meta”). The branching ratios 
are obtained from Tables lA~4l and lA.5l This reaction produces the 
Fi, Ai, and E species of ND^ in ratios 9:9:2. 


Fig.A.l. Branching ratios of the reaction - 1 - NH 3 —> 
(NHg)* —> NH 4 + H 2 when the reactants are of the Ai (“ortho”) 
and E (“para”) symmetries. 


1/1 


Ai 0 A 

A; 0 A 

Ai 0S 

F] 0A 

E®B 

F] 0A 
Ai 0 A 

E®A 

El 0A 
Fi 0S 

Fig.A.2. Branching ratios of the reaction + ND 3 —> 
(NDg )* —> ND 4 -H D 2 when both reactants are of the Ai (“meta”) 
symmetry. 
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Table A.4. Statistical nuclear spin branching ratios of the reaction Dg ^ D 3 -H D 3 . 



Table A.5. Statistical nuclear spin branching ratios of the reaction Dg ^ D 4 -H D 2 . 
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